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SECTION I 


INTRODUCTION 

Atomic polarizabilities describe the changes in the charge distri 
bution of an atom when it interacts with an electric field and because 
of its close relationship to dielectric constants, the study of dipole 
polarizabilities has a lengthy history. Interest in atomic shielding 
factors is of more recent origin and was stimulated by attempts to 
determine nuclear quadrupole moments from observations of hyperfme 
structure. To a first approximation, atomic hyperfine structure is due 
to the interaction between the nucleus and the valence electrons and 
the inner closed shells of electrons act merely to screen the nuclear 
charge and so to modify the electrostatic field in which the valence 
electrons move. However, another type of shielding occurs because the 
valence electrons distort the spherical symmetry of the closed shells 
and the distorted shells then interact with the nucleus. Alternatively, 
the nuclear multipole moments induce equivalent moments in the charge 
distribution of the closed shells and these moments then interact with 

the valence electrons . 
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Although they refer to different phenomena, polarizabilities and 
shielding factors are examined together in this review since the math- 
ematical procedures for calculating them are very similar and a larger 
body of experimental data becomes available for assessing the accuracy 
of the methods used. 

1.1 DEFINITIONS 

When an atom is placed in the field of an external charge Z' , it 
is polarized and the resulting distribution of charge can be character- 
ized by a series of induced electric multipole moments, each of which 
is proportional to Z' provided Z' is small. The induced dipole moment 
is related to the electric field of the charge through the atomic dipole 
polarizability, a d . The induced quadrupole moment, a q , is related to 
the gradient of the electric field of the charge through the atomic 
(field-gradient) polarizability * Similar relations involving higher 

order derivatives of the electric field apply to the higher pole polar- 
* k 

izabilities . 

The electric field of the external charge is also modified by the 
charge distribution of the atom and the dipole shielding factor, (3^, 
can be defined as the ratio of the change in the electric field at the 
nucleus due to the charge distribution of the atom to the electric field 
at the nucleus due to the external charge alone. The quadrupole shield- 
ing factor, 7 , is the ratio of the change in the gradient of the elec- 
trie field at the nucleus due to the charge distribution of the atom 

Occasionally a different meaning is attached to quadrupole (and higher 
pole polarizabilities). The uniform part of the electric field induces 
by a second order dipole-dipole coupling, a quadrupole moment which is 
proportional to Z /2 . The constant of proportionality is also referred 
to as a quadrupole polarizability. 
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to the gradient due to the external charge alone. Similar definitions 
involving higher order derivatives of the electric field apply to higher 
pole shielding factors. Shielding factors can be defined alternatively 
as the ratio of the multipole moment induced in the electron charge* 
distribution to the nuclear multipole moment which gives rise to it. 

When the atomic system is not spherically symmetric, the distor- 
tion of its electron distribution by the electric charge will depend 
upon the orientation of the atomic system. In such cases, it is con- 
venient to average the polarizabilities and shielding factors over all 
orientations . 


1 . 2 QUANTAL FORMULAE 


Suppose the unperturbed atom has N electrons with position vectors 

r referred to the nucleus as origin and let Z be the nuclear charge . 

~i 

The unperturbed Hamiltonian is given in atomic units by 


N 


H = - 


I 


i=l 



( 1 ) 


and the unperturbed eigenfunction \|r (r) , when r denotes collectively 

all the position vectors r., satisfies the Schrodinger equation 

~1 

(H ■ V *o ( u = °* (2) 

E being the unperturbed eigenvalue, 
o 
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Suppose also that the external charge Z' is located at x' and that 


x' is large. Then the interaction potential 


V(r, £') 




(3) 


may be expanded in the form 


N oo k 

v( £ , x') - Z ' K \' N) - z ' Z I V cos V 

i~=i if=i r 


(4) 


where (r., 0 i » (t> i ) are the spherical polar co-ordinates of r^ the polar 

axis being chosen such that the co-ordinates of x' are (r 7 , o, o) and 

P (cos 9 ) is the Legendre polynominal of order k. The constant 
k i 

spherically symmetric part of V(r, r 7 ) is of no interest here and we 
replace (4) by 


V(r,r') 



P, (cos 9 . ) . 
k i 


(5) 


Let \(r (r | r 7 ) be the wave function of the system when a charge 
Z' is placed at r 7 . Provided the unperturbed state is not degenerate 
(as we shall assume), the perturbed wave function may be written in the 

form 00 (k) 

i(f (r | r 7 ) = >|r q (r) + Z' ^ + \ + 0(Z 7 ), (6) 

k=l 


t 


(k) 


being chosen so that 




(k) 


V = ° 


(7) 
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for all values of k. The requirement (7) ensures that f is normalized 
to the order of Z\ 


We shall define the electric multipole moment of order 2 of a 
charge distribution as 


a 





(cos 


\ 

9 .) t • 
1 / 


( 8 ) 


To the order of Z' , (8) is 


fl 2 L 



1 

,k + 1 

r ' 



The electric field strength at the nucleus due to the charge Z' is 
o 

-Z / /r / and the dipole polarizability is given by 


a 


d 


= 2 




r . 
l 


(cos 0 .) \lr 
x i o 


( 10 ) 


3 

The gradient of the electric field at the nucleus is -2Z / /r / and the 
quadrupole polarizability is given by 


a 


q 


= 2 


N 



i=l 



P 2 (cos 0 ± ) 



(ID 


the factor of 2 being inserted so that (11) conforms with the defini- 
tion of employed in most of the literature (c.f. Sternheimer 1954). 

We are assuming that the wave functions are real. If they are not, 
the real part of the expression is to be taken. 
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In general, the 2- pole polarizability is given by 

N 


a 2 L = 2 


( 1 : 

\ i=l 


J P L (cos ep i o 


( 12 ) 


in units of a 2L+1 , as being the Bohr radius 5.292 x 10 9 cm. With this 
o 

definition, the polarizabilities of any atom in its ground state are all 
positive quantities . 

The electrostatic potential at the nucleus due to the atomic 
charge distribution and to the external charge is 

N 


$ (r 


= ■ X ^ + f 7 


(13) 


L=1 


and the expectation value of the electric field at the nucleus is 
accordingly 

P (cos 0.) 


( N P^ (cos 0. 

+• y — 


which to first order in Z is 


(14) 


2Z / 


L r7 *r ( I ^ 

k=i \ i 




i' 


(15) 


The dipole shielding factor is therefore 


j3oo = 2 


( m V p i (cos V ■ 

( +i -1 *o 

V i=l i 


(16) 
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Feynman (1939) has shown that the force on a nucleus of charge Z 

is Z times the electric field at the nucleus due to all the electrons 

plus the fields from the external charge. The force of an ion charge 

2 

Z' acting on an ion of net charge (Z - N) is simply -(Z-N) Z'/r' in 
the positive r' direction and the corresponding total electric field is 
-(Z - N) Zr / ^. Comparing with (15), it follows that 


fi _ s 

p " z ’ 


(17) 


an argument first published by Sternheimer (1954). The dipole shield- 
ing factor is consequently a known quantity. Its interest lies in the 
fact that (16) contains the first-order perturbed function ^ (r) which 
is required for the evaluation of the dipole polarizability. The rela- 
tionship (17) provides an assessment of the accuracy of the calculated 
first-order functions. 


The expectation value of the gradient of the electric field at the 
nucleus is N 




2P 2 (cos Q /) ^ 2 z ' 

T|r 


i=l r i 


r /3 


(18) 


which to first order in Z / is 


N 


2Z ' 


k + 1 


f y 

i ’ / , 


k=l 


2P 2 (cos 0.) 


i=l 


r . 

l 


2Z' 


,3 


(19) 


7 




which is dimensionless. With this definition, a positive value of 
corresponds to shielding and a negative value to antishielding . 
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SECTION II 


PERTURBATION THEORY 
2.1 THE FIRST ORDER EQUATION 

(k) 

The description of the determination of (r) can be simpli- 

fied by treating each term of V (r, irO 3S a separate perturbing poten- 
tial. Thus, writing 

N 

v fc (r) = - Yj r i kp k (cOS 0 i )j (22) 

i=l 

the perturbed wave function ^(r) can be expanded in the form 

\(£> = t Q (r) + Z' 1^°° (r) + Z' 2 * 2 (k) (r) + 0(Z /3 ) (23) 

and the perturbed eigenvalue in the form 

= E + Z' e. (k) + z /2 f 9 (k) + 0(Z /3 ). (24) 

k o 1 ^ 

Then the wave equations of the perturbed system 

(H + Z^ k - & k ) t k = 0 (25) 
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may be replaced by the sequence 


(H - E ) i =0 
o o 


( 26 ) 


00 


00 


(H - E o ) V + (v k - V ' = ° 


(27) 


(H - E ) * (k) + (v k - ^ (k) " <5 ' 2 (k> = 0 


(28) 


From (27) 


6 l (k) - < V v k V 


(29) 


and from (28), 

* (k) 


( v k * o ) 


(30) 


But from (12) 


(L) 


a 2 L - -2 ( Y v L + o ) 


(31) 


so that 


f (L) 
6 2 


" 2 a 2 L 


(32) 


The shielding factor (21) may be written 


7 2 L = 2 ( 'I' 1 (L) , v' L * o ) 


(33) 
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where 


v 


/ 

k 


(r) 



P, (cos 0 . ) 
k i_ 

k + 1 


(34) 


It follows from symmetry considerations that (33) is identical to the 
alternative form 


7 2 L = 2(+ 1 (L)/ , v L t o ) 


(35) 


where 


(H - E q ) + (v' k - = 0 


(k) 


(36) 


such that 


V -0 < 37) 

and 

- ( * 0 . v' k * o ) (38) 

A formal proof of the identity of (33) and (35) follows from the 
use of (26) to yield 


( *■ 


(L) 


V = 


( 


(L) ' 


H - E 


V. 


(L) 


) (39) 
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and (36) to yield 


-( t 1 (L) | H - E o | ^ (L) ) = ( v' L Hr o ) 


(40) 


2.2 THE HYDROGEN ATOM 

As a demonstration of the theory of § 2.1, we consider a hydrogen 

(k) 

atom in its ground state. Equation (26) for ^ can be simplified 
by writing 


^ (k) (K) = (O (O • 


(k) 


(41) 


(\r\ 

Then J (r) satisfies the equation 


N 


N 


V 

r(k) i o 


(k) 


X v i 2 £ (k) - £ v i £W • + \-V''-° 


i=l 


i-1 


For a hydrogenic atom in its ground state 


(42) 


t (r) 

o ~ 


= (Z / jt ) 2 exp (-Zr) 


(43) 


and (42) reduces to 


_ I V 2 f( k ) + z 7 ^ - r k P (cos 0 ) = 0, 

2 dr k 


( 44 ) 


with solution (Dalgarno and Lewis 1955) 
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P (cos 6) . 

K 


(45) 


f 



k + 1 
r_ 

k + 1 


r 

kZ 


From (30) the multipole polarizability is 




(2L + 2} (L + 2) 

2 2L +1 Z 2L + 2 l(l + 1} 


(46) 


and from (33) the multipole shielding factor is 


7 2 L L(L + 1) Z 


(47) 


The solution (45) for k = 1 was first obtained by Waller (1926), 
who calculated the dipole polarizability of atomic hydrogen. 

2.3 UPPER BOUNDS 

Except for one-electron systems it is not possible to solve (27) 
exactly. It may be solved formally by substituting the expansion 

/ 

+i (k) <t> - S a t <k> *t w (48) 

t 

where \|r (r) is the eigenfunction of the t th excited state of the 
t ~ 

unperturbed atomic system, satisfying the wave equation 


(H - E t ) * t (r) = 0 


(49) 
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and the prime on the summation ir.dicatas that the term t = 0 is excluded, 

(k) 

Then the expansion coefficients a^ are given by 


00 


( v v k v 

E t - E o 


(50) 


and 


a 2 L = 2 


/ hy v l Vj 2 

E t - E o 


(51) 


which is clearly positive for all ground state atomic systems. 

The convergence of the summation over t as the number of terms is 
increased is usually slow, as Tillieu and Guy (1954) have explicitly 
demonstrated for atomic hydrogen. If transitions to continuum states 
are unimportant, the convergence can be made more rapid by using the 
equivalent formula (Lennard - Jones 1930) 


a 2 L = a 2 L + 2 


(e i ( l) - y i ( v v l v 1 2 

(Ei (L) - E o ) (E t 


V 


(52) 


where 


^ < E i a> - E o> 


2 2) 

( V V L V ' ( V V L *o ) J’ 


(53) 
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and the prime on the summation :*.n*...'_.r;at33 that the term t 0 is 

(k) 

Then the expansion coefficients a t are given by 


excluded. 


a 


t 


00 _ 



v k +o> 


- E 


(50) 


and 


= 2 


5 Ijv v l v 


(51) 


which is clearly positive for all ground state atomic systems. 


The convergence of the summation over t as the number of terms is 
increased is usually slow, as Tillieu and Guy (1954) have explicitly 
demonstrated for atomic hydrogen. If transitions to continuum states 
are unimportant, the convergence can be made more rapid by using the 
equivalent formula (Lennard - Jones 1930) 




u ^ 

C^L +2 

t 


( e i ( l) - v i (y v L vj 2 

(Ei (L) - V < E t - V 


(52) 


where 


= (e/ L >-E) 
1 o 


( V V L 2 V ‘ ( V V L 


(53) 
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E being the eigenvalue of the lowest excited state for which 

v W ) does not vanish. With few exceptions, neither (51) nor (52) 

L r o 

provides a useful means of computing ^ is valuable to note 

that, since - E t is negative is an upper bound 


a 2L - a 2L 


(54) 


2.4 OSCILLATOR STRENGTH FORMULA 


Formula (51) may be used to derive accurate values of the dipole 
polarizabilities for the few atomic systems for which reliable oscillator 
strengths are available. The electric dipole oscillator strength corres- 
ponding to the transition from the state i|r Q to the state is defined 
by 


f 


ot 



- V 


2 

<V E V 

i 


(55) 


and in terms of it, 

/ 

“ot - S £ ot /(E t - V 2 ' <56) 

'V' 

t 

Dalgarno and Kingston (1959) have used (56) to predict the polari- 
zabilities of the alkali metals and their results are given in Table 1. 
The error is unlikely to exceed 10%. Their predictions disagree with 
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early measurements but have been confirmed by the recent experiments of 
Salop, Pollack and Bederson (1961) whose results are also included in 
Table 1. 

1 3 

The dipole polarizabilities of metastable helium in the 2 S and 2 S 
states (Dalgarno and Kingston 1958) and of the negative hydrogen ion 
(Bates and Lewis 1955, Geltman 1962, Dalgarno and Ewart 1962) have also 
been obtained by summing the individual terms of (56) but because of the 
contribution from transitions to the continuum, the method is rarely 
convenient. Donath (1962) has essentially used (56) to predict the 
dipole polarizabilities of F _ , Ne and Na + , the wave functions of a num- 
ber of discrete states being determined by variational methods. The 
results are not encouraging and a direct variational attack on the com- 
plete summation is usually to be preferred. 

2.5 APPROXIMATE FORMULAE 

Approximate evaluations of (51) can be made by replacing E t - E q 
in the denominator by a mean excitation energy E, yielding 

°2L ~ I {«V V L 2 *o> - <V V L *o> 2 } • (5?) 

A common choice for E is the ionization potential of the atom and 
the resulting formula usually yields values of w hi- c h are correct to 

within an order of magnitude. 
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TABLE 1 


-24 3 

DIPOLE POLARIZABILITIES OF THE ALKALI METALS IN UNITS OF 10 cm 



Li 

Na 


K 

Rh 

Cs 


Formula (56) 

24.4 + 2.4 

24.6 + 2.5 

41.6 

+ 4.2 

43.9 + 4.4 

53.8 + 

5.4 

Measured 

20 + 3.0 

20+2.5 

40 + 

5.0 

40+5.0 

52. 5 + 

6.5 


A more interesting choice of E is such that 

O / 2 _ 2 

O (E t - E o ) I (t t , v L i|f o ) I = E O |(i t , v L i o )| (58) 
t t 

(Vinti 1932). We shall demonstrate in^3 that it leads to a lower bound 
for a 2L- 

It is apparent from the form of (57) that (56) is valuable for 
those systems in which 

( V V L V 2 >;> V L V 2 >:> V L ^o) 2 ' ‘ * (59) 

(Vinti 1932), as is often the case. 

The procedure of replacing E^_ - E q by a mean energy E also yields 
approximate formulae for shielding factors 
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but they may be very misleading; the condition analogous to (59) is not 
satisfied, transitions to continuum states being of great importance. 
This is reflected in variational calculations by the sensitivity of 
shielding factors to the assumed trial functions. 
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SECTION III 


VARIATIONAL METHODS 


3.1 APPROXIMATE FORMULAE 

As (32) demonstrates, electric polarizabilities are directly related 
to the change in energy due to the perturbation and they may be deter- 
mined therefore by application of the conventional Rayleigh-Ritz varia- 
tional methods. Many of the early studies proceeded by substituting a 
trial wave function \|r for the perturbed system into the expression for 
the total energy 


£= (t t |H + v| t t ) 


(61) 


and minimizing, the final expression being expanded in powers of the pre- 
turbing potential. Xt may be shown (Slater and Kirkwood 1932, Dalgarno 
and Lewis 1956) that provided is the exact unperturbed eigenfunction 
this procedure is equivalent to the simpler one of minimizing the func- 
tional 


& a) = 

2 



| (H - E | (L) ) + 2 (t (L) | 

t t 


V L -V l) l +„> (62 > 
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with respect to a trial first order perturbed function • In prac- 

tice \lf is not the exact eigenfunction and the former method introduces 
o 

spurious terms, which if retained may seriously affect the derived polar- 
izability. If they are ignored, the procedure is again equivalent to 
the use of (62). 

The variational trial wave function 


*1 


t 


(L) 


00 = 



V£> 

O ~ 


(63) 


if of special interest. Writing 


V L = V L " e i 


(L) 


(64) 


it follows that 


\ - (V L * o |H 


-<*o’ V L V 


VX V 


(65) 


and 


£ CD = 
2 


■ ( V V L 2 

(v L * 0 |H - EJ v L * o ) 


(66) 


Thus 


a 


2 <V \ v 2 


2L (v_ 1 

L o 


H 


V L V 


( 67 ) 


is a lower bound for Q^t • 
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It is a simple matter to show that the denominator of (68) is equal 


to the infinite summation (59) 


S (E t - E o )| (V v L i o )| 2 = (v L ^ |h - E o |v l * o ) 


( 68 ) 


(Vinti 1932 , Dalgarno and Lewis 1957), thereby proving that the choice 
of E described by (59) leads to a lower bound. 

The denominator of (68) can be written in the simple form 


N 


*„I H - E ol v L *o> - - <V L V Z V i V L ' V i V 


(69) 


i=l 


For L = 1, integration by parts shows that 


(Vi * o |H - E o |v 1 * G ) = \ N 


(70) 


(this being the Thomas-Kuhn oscillator strength sum rule when taken 
with (68) ), so that 


a d> 


4 5 V 1 
1 

N 


(71) 


(Kirkwood 1931, Vinti 1932). For heavy atoms, the lower bound given by 
(71) j^s very much less than the actual value and a better indication is 
of the possible magnitude of is provided by replacing N by the number 
of electrons in the outmost shell. 
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For most atomic system in states of zero orbital angular momentum 


(but not all), 



I E 


r . ° r , 

'V J 


♦ ) < o 

o 


and (70) can be replaced by the expression 


N 


a . 


> ±_ (f y r . 2 ¥ )‘ 

- 9N V o’ l_, 1 o / 


i=l 


(72) 


(73) 


The diamagnetic susceptibility of a gas is 


X = 


2 2 
e L a 

o 


6 me 



(74) 


where L is Loschmidt s number (c.f. Van Vleck, 1932) and Kirkwood (1931) 


proposed the semi-empirical formula 


he 2 a 1/2 ,,, 

X 1 (Na) 1/2 , 

4 me 


(75) 


later modified by Vinti (1932) to take account of the r , . r . terms 
and by Buckingham (1937) to give more precise significance to the value 
of N. The formula is surprisingly successful (Brindley 1933 5 Buckingham 
1937). 
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Simple formulae can also be derived for higher order polarizabilities. 


Thus for atoms in S states, 



(76) 


(Dalgarno and Lewis 1957) so that 


a > 

q - 



(77) 


Dalgarno and Lewis (see also Stewart 1961) have suggested several formulae 
for estimating the order of magnitude of which are based on (77). 

They suggest, in particular, 


a. ~ 

q 




(78) 


which has the merit that it does not require a knowledge of any wave 
functions . 

The reliability of these various approximate formulae depends upon 
the degree with which condition (59) is satisfied and there is little 
value in carrying through a similar development for shielding factors. 
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An interesting generalization of (63) has been proposed by Pople 
and Schofield (1957). They use as trial function for the calculation 
of dipole polarizabilities 

N 

t 1 (L) = U) f(r i } C ° S 9 i (?9) 

i=l 

where f(r.) is allowed to vary arbitrarily. They show that if ^ Q (£) is 
represented by the Hartree-Fock approximation the optimum form of f(r^) 
is the solution of a simple differential equation which may be integrated 
numerically. The computation involved is less than that of the Stern- 
heimer procedure (Para. 4.4) and the method does not encounter any diffi 
culties arising from exchange effects, (79) being a properly antisym- 
metrized function. For argon, Pople and Schofield obtain a dipole polar- 

izability of 2.03A 3 compared to 2.40A 3 given by the Sternheimer proced- 

°3 

ure (Kaneko 1959) and to 1.64A given by extrapolation of refractive 
index data (Dalgarno and Kingston 1960). Because of its simplicity, the 
method merits further application though it is less accurate than the 
uncoupled and coupled Hartree-Fock approximations (Para. 4). 

3.2 TWO-ELECTRON SYSTEMS 

More complex trial wave functions have been used for the two-electron 
systems H _ , He and Li + . It is clear from the work of Hasse' (1930, 1931) 
and of Barber and Hasse' (1937) that a d is sensitive to the adopted repre- 
sentation of the unperturbed wave function \|f Q but that provided the form 
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of the trial perturbed wave function is carefully chosen, conver- 

gence of is obtained with a small number of parameters. Using tho 
six-parameter representation of due to Hylleraas (1929) and the trial 
wave function 

(r L , r = ^( ar L + hr i 2 + CT] 3 ) cos> 0 i + ( ar 2 + br 2 + Cr 2 ^ 

cos ©2 ^ (^0) 

Baber and Hasse 1 obtained (0,2045 + 0,0005) x 10 cm for the dipole 

polarizability of helium. More recently Schwartz (1961) has used an 18- 

parameter representation of and a much more flexible trial function 

2^ 3 

and he computes a value of (0.2050 + 0.0001) x 10 cm . This 

value does not agree with the most recent directly measured value of 

rs / n 

(0.2068 + 0.0002) x 10" cin (Johnson, Oudemans and Cole 1960) but it 

-24 3 

is in harmony with the value of 0.2051 x 10 cm obtained by Dalgarno 
and Kingston (1960) from a careful extrapolation of refractive index 
data . 

A different type of trial function has been investigated for helium 

^ , 1 

by Abbott and Bolton (1954), which is essentially a sum of excited P 
state wave functions. This choice must share some of the convergence 
difficulties associated with the evaluation of (56) and compares unfavor- 
ably with (80) . 
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Baber and Hasse' (1937) have also obtained a value of 0.02862 x 
10~ 24 cm 3 for the dipole polarizability of Li + which is again based on 
the 6-parameter representation \|r and a trial function (80). The value 
obtained from an analysis of the spectral term defects (Mayer and Mayer 

o / n 

1933, Sternheimer 1954) is 0.0235 x 10" cm but the variational ly 
determined value is undoubtedly the more accurate. 

Schwartz (1961) has also calculated the dipole polarizability of H 

using wave functions similar to those he adopted for helium. He finds 

that a = 31.4 x 10~ 24 cm 3 * which agrees well with that obtained, CC^ = 
d 

30.2 x lO - ^ cm^ (Geltman 1962, Dalgarno and Ewart 1962) by evaluating 
(56) from the photo-detachment cross sections computed by Geltman (1962). 

The only refined calculation of a quadrupole polarizability and a 
quadrupole shielding factor appears to be that of Dalgarno, Davison and 
Stewart (1960). They employ the Hylleraas representation of t Q for helium 
and a trial function similar to (80) but with an additional parameter 
and obtain a = 0.0942 x 10 _4 ° cm 3 and 7oo = 0.397. As they demonstrate, 
a is rather sensitive to the adopted form of \|r and it is probable that 

q 

their value is too small, the uncertainty being about 5%. The value of 
7 is probably accurate to within an error of 24. 

Detailed variational calculations of the type applied to two-electron 
systems are scarcely feasible for more complex systems and for them it is 
natural to seek an extension of the Hartree-Fock scheme for unperturbed atoms 

— = — — — ,t 2 A- 3 

*The value published by Schwartz is 26.8 x 10" cm . I am informed by 
Dr. L. M. Branscomb that this value is incorrect and that Schwart s 
calculations actually yield the value quoted in the text. 
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SECTION 4 


THE HARTREE-FOCK APPROXIMATION 


An extension of the Hartree-Fock approximation scheme to the deter- 
mination of polarizabilities and shielding factors has been carried out 
by Dalgarno (1959), Kaneko (1959), and Allen (I960). The essential 
features are contained in earlier discussions by Peng (1941) and by 
Temkin (1957). Because of its close connection with earlier work, we 
begin with a description of the Hartree approximation in which electron 
exchange effects are ignored. 


4.1 THE HARTREE APPROXIMATION 


In the Hartree approximation, it is assumed that the unperturbed 
functions may be written in the form 




N 


n 



(O 

~1 


(81) 


where the orbitals u. ^ (r.) are eigenfunctions of the equations 


t u i (0) ( ^i ) = ( 1/2 V i 2 + W) U i (0> (£ i } = e i U i (0) 


H 


(82) 
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the self-consistent potential V (r^) being defined by 


(r k ) I 2 


Z V f Uk — 

W " r ‘ + £ ■* l£i • £■ 


dr k (83) 


and the summation being over the occupied orbitals. 


The orbitals u. (r.) are written as the product 
1 


, _ P (r, £. \ r. ) 

u (o) (r.) = - - — - — — Y (6 ,4>.) 

1 i ii 


(84) 


where n , £ and m are respectively the principal, azimuthal and magne- 
i i l 

tic quantum numbers specifying the orbital, P q (n^i^|r^) is the radial 
wave function normalized so that 


CO 

/ WJV W f il r i ) dr = Bn iV 


(85) 


an< j y (0 ) is the spherical harmonic of order £ . Then (82) 

£ ,m. i i x 

i i 

reduces to 



+ V. (r. ) + 


l i 


£.(£. + 1 ) 

i i 

2r . 2 



P (n,£. r, ) = 0 
o v 1 1 1 1 


( 86 ) 


28 



The total unperturbed function \|/ Q (r) is a solution of 


(H 7 - E o ') t OO = 0 


( 87 ) 


where 


H 7 




and 



(89) 


The first order perturbed wave function is written in the form 


K u. (l> (r,) 


^ * e<s> 

1=1 X 


(90) 


where we have suppressed the superscript L. Corresponding to (90), the 
electric multipole polarizability is given by 
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t 


a 


2L 


= 2 



(91) 


and the multipole shielding factor is given by 


7 2L 



p (cos 6 . ) 
L 3- 

L + 1 
r . 



(92) 


The procedure of Pople and Schofield (1957) is equivalent to (90) with 
the restriction that u ^/u^^ be independent of the electron shell. 

4.2 THE UNCOUPLED HARTREE APPROXIMATION 

There are two procedures for specifying the perturbed orbitals u^ 

(r.) (Dalgarno , 1959). In the simpler procedure, the radial parts of 
~i ° 

u and u ^ are allowed to undergo independent but otherwise arbit- 

i i 

ra ry variations. This is equivalent to the assumption that 


(H / - E q / ) >|/ L (r) + v L t Q (r) = 0. 


(93) 


Substituting (90), it follows that 


(H i ' £ i (0)) U i (1) ( -i } + V L ( ^i ) U i W ' ( ~i } “ ° 


(o) 


(94) 
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where 


(r\) = r i L P L (cos 6 1 ) - (^ u i ( ’°' ) > r i Lp L ( cos 0 i ) u i^y » (95) 


it being convenient to require that 


^ a) . ^ ■ ° 


(96) 


4.2.1. Expressions for an< ^ ^2L* Foley, Sternheimer and 

Tycko (1954) pointed out that (94) can be reduced to a finite set of 
uncoupled radial equations. Thus, generalizing their procedures, we 
express the inhomogeneous angular part of (94) as the finite sum 


P. (cos Q . )Y (0 , ) 

1 v i Jb .m, l i 

i l 



2 £. + 1 \ ly ^ 2 

TT7r-riJ c(L Vi -om.) 


C(L ;oo) 
'll 



( 0 i> 


>i> 


(97) 


where C(i is the Clebsch-Gordon coefficient in the notation 

of Rose (1957). Accordingly, we may write 


u . 

l 


( 1 ) 


(r.) 

~i 



2 h + 1 \ 

2x ' + i y 


1/2 


C (LlL-gy ; 0m i ) C(Li j .i i / ;oo) 
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( 98 ) 


V n iW |r 


x / m i ( 6 i’ V r 


and ( 94 ) becomes 


r h 2 £ ± ( V + 1} 

1/2 A + V 4 (r 4 ) + 


dr 


2 i' i 


2 r 


i <o) } WuV 


00 N 

f rL - 5 f.i / P 0 (n t 4 i |r) rl P o <Vil r > dr ) Wj 0 
V 1 o 


for all values of £^ included in the sum ( 97 ). In solving ( 99 ), 
must be taken to insure that f P 1 (n i i i |r) P Q (n i i i |r)dr = 0 . 


The multipole polarizabilities become 


C£„ =2 

2 L 


r — , / 2H , . + 1 X 2 

Y y l — \ C(Li J.^oo) x ynJ.jf.O 

n.£ £.' \ 2i i + 1 / 


z 

m. 


C (LM^i 7 ; 


om. ) 

l 


2 


= 0 
( 99 ) 


care 

( 100 ) 


( 101 ) 
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where 


Wi 


>v> 


J P l (n i £ i ;£ i / \r) r L P o (n.£.lr)dr 

o 


( 102 ) 


and the summation is over the occupied orbitals. Averaging over all 
orientations of the undisturbed atomic system, (101) simplifies to 


a 

2L 


Z I 


a(n i £ i ) 

(2L + 1) 


C(L £.£.' ;oo) 

XI 




(103) 


n. £ . 

l l 


Jo.' 

i 


where a(n £ ) is the number of electrons in the (n £ ) shell. The 
i i 

parity C-coef ficient ;oo) vanishes unless £ ± ' lies between L 4- £ 

and |L - £ I and unless L + £ . + £ ' is even. Its value is 

1 l 1 IX 


M 2 


, . . . ,1/2 (L + £ . 

C(L£ i £ i / ;oo) = (-1) i 


v> 


2 £ . ' + 1 
x 


L + £ . + £. / + 1, 

X X 


X 


(104) 


t (L + £ . + £.') 

V X X 

(L + £ , - £.') t (L - £. + i/) T (" L + £ i + 

N X X 1 1 ■ L 


where 


t(x) = (l/2x) ! / ^x ! 


(105) 


(Racah 1942) . 
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Straightforward evaluation of (104) yields for closed shell systems, 


the dipole polarizability formula 


a d 


4 

3 


V < (f + 1) I, (ni;i + 1) + il 1 (n £;£ - 
n£ ^ 



(106) 


the quadrupole polarizability formula 


a 

q 


6 

5 


z 

n£ 


(J, + 1) (£ + 2) 
(2^ + 3) 


I 2 (ni;i + 2) + 


21(1 + 1) (2i + 1) 
3 (2i - 1) (2£ + 3) 


I 2 (ni;i) + 


(i - l)i 

(2i - 1) 


I 2 (n^ 



(107) 


the octupole polarizability formula 



n£ 




(i + i) ce + 2) (i + 3) 


(2i + 3) (2£ + 5) 


I 3 (ni;i + 3) + 


3l£l +1) (i+ 2) 

5 (2i - 1) {21 + 5) 


I_ (a 6 ;£ + 


1 ) 


(i - 1) i(i + 1) 
(2£ - 3) (2£ + 3) 


I 3 (ni;i - 1) 


(J - 2 ) g - 1 ) 
(2i - 3) (2i - 1) 


I 3 (nj? ;£ 



(108) 


and the hexadecapole polarizability formula 
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a 


12 



a +i) a + 2 ) (£ + 3 ) (i + 4) 

(24 + 3) (24 + 5) (24 + 7) 


I 4 (ni;i + 4) 


4 4(4 +D (i + 2) (4 + 3) 

7 (24 - 1) (24 + 3) (24 + 7) 


I 4 (n4;4 + 2) 


18 (4 - 1) 4(4 + 1) (£ + 2) (21 + 1)__ 

35 (2i - 3) (24 - 1) (2£ + 3) (2i + 5) 


I^(n £ ;£) 


4 a - 2) (i - 1) £(i + 1) 

7 (2J - 5) (2£ - 1) (2£ + 3) 


I 4 (ni;i - 2) 


(£ - 3) {£ - 2) (i - 1) _ 
(2i - 5) (2i - 3) (2£ - 1) 


I 4 (ni ;£ 


- 4 >) 


(109) 


Some of the terms of (106), (107), and (109) for small values of £ 
have been derived previously by Foley, Sternheimer and Tycko (1954) and 
Sternheimer (1954, 1961). 

The multipole shielding factors are given by expressions identical 
to (101) - (109) except that Kniji') is replaced by 

00 

J(n £-,£') = / P 1 '(n£;i'|r) L ^ ~ P 0 (ni|r)dr. (110 
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Alternatively, I(ni;i / ) may be replaced by 


J'(n &•,£') 


I 


P 1 / (n£;£ / | r)r L P Q (n^|r)dr 


(111) 


where V ^ (ni, \r) is the solution of 


{■ 



+ V.(r) 


£'(£' + 1) 

+ 2 
r 



{n£-,£' 


r) 


+ 


{a 1 




0° -v 

J p „M r) ;t L TT p o (nf l r ) dr j 


P q (ni|r) = 0. 


( 112 ) 


4.2.2. Solution of First Order Equation. Equation (99) may be 
solved by direct numerical integration and Sternheimer (1954, 1957) has 
described in detail a suitable computational procedure. He points out 
that the labor involved in a direct evaluation of V^(r^) can be avoided 
by using the relationship 


V. (r. ) - e. 
i ~i i 


(o) _ 


= 1/2 


P(n.i. I r. ) 

l i* i 


d 2 P(n i ^ i [r i ) 

2 2 f 

dr. r. J 


(113) 
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The method of numerical integration usually involves several iterations 
and is tedious to apply. 


Equation (99) may also be solved by variational methods. It is 
instructive to return to Equation (94) and to construct the functional 





+ 2 


/ 



V L ( Ii > 



(114) 


Minimizing 



with respect to 


the trial function 


u . 

l 


(i) 


<£i> 




u.<°> (r.) 
1 


(115) 


leads to the dipole polarizability formula 


a 


d 



(116) 


which was derived by Heilman (1935) and by Buckingham (1937) by super- 
ficially different methods. Several of the other more complicated form- 
ulae derived by Buckingham (1937) and by Bravin (1953 5 1954) also follow 
directly from (114) using the more flexible trial functions 


u . 

i 


( 1 ) 


(r.) 

~1 


(a + br. + cr 2 ) v (r ) u ^(r ) 
111 ~1 1 


(117) 
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It is clear from (98) that, except when i - 0, (117) has the wrong 
angular dependence. 

The simplest trial orbital with the correct angular dependence is 
given by (98) with 

P 1 (n i ^ i ;i i / |r) = X ( n jA r p 0 ( 118 > 

but it leads to negative polarizabilities (except when & = 0). With 
more complicated trial wave functions, it is possible to obtain solutions 
of (94) of comparable accuracy to those obtained from numerical integra- 
tion and with less labor. It has been common practice in recent years 
to select trial functions of the form 


Wi ! Vl r > - Y. 


a S r S P (n.i . | r) , 

Oil 1 


(119) 


a choice which has the advantage that by using (84) the potential V^(r) 
can be removed from the calculation. It has the disadvantage that it 
forces nodes into P L (n i | r) which should not occur (Sternheimer 1959). 
The effect on the derived polarizabilities is usually not serious, partly 
because the polarizabilities are determined by the distortion at large 
distances from the nucleus, but many of the variationally determined 
shielding factors are of uncertain accuracy. 
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4.2.3. The Helium Isoelectronic Sequence. The only systems to 


which the Hartree approximation may properly be applied are those com- 
prising the helium isoelectronic sequence. The dipole and quadrupole 
polarizabilities and shielding factors have been obtained by numerical 
integration of (99) (Sternheimer 1957, 1959) and by variational methods 
using (117) (Das and Bersohn 1956, Wikner and Das 1957). There are 
small discrepancies between the two sets of values which can be attri- 
buted to the different representations employed for the unperturbed 
orbitals. The results are given in Table 2 which includes also the most 
accurate estimates available of the various quantities. 

With the notable exception of H~, the accuracy of the uncoupled 
approximation is good. It improves rapidly with increasing nuclear 
charge and it is better for polarizabilities than for shielding factors, 
suggesting that the approximation provides a more accurate description 
of the distortion in the electron density in regions far from the nucleus 
than in regions near to the nucleus. 

For more complex systems, the influence of exchange effects requires 
cons id er at ion. 

4.3 EXCHANGE EFFECTS 

Buckingham (1937) found that (116) yielded dipole polarizabilities 
in better agreement with experiment if the unperturbed Hartree wave 
functions were replaced by the unperturbed Hartree-Fock wave functions 
and following Sternheimer (1954) it has been common practice to similarly 
modify the uncoupled Hartree approximation. Thus, in (99) and in (113), 
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X 


TABLE 2 

DIPOLE AND QUADRUPOLE POLARIZABILITIES AND SHIELDING FACTORS OF 
THE HELIUM ISOELECTRONIC SEQUENCE 

a d (10' 24 cm 3 ) a q( n " 40 Cm5) 


Uncoupled Coupled Uncoupled Coupled 

Approximation Approximation Accurate Approximation Approximation Accurate 


h“ 

13.4 

-11.8 

31.4 

66.5 

- 

- 

He 

0.220 

0.196 

0.205 

0.0979 

0.0965 


Li + 

0.0304 

0.0281 

0.0286 

0.00472 

0.00464 

- 

Be’ H ” 

0.00820 

- 

- 

0.000637 

- 

- 



Poo 


7 

'oo 



H _ 

3.25 

2.00 

2.00 

- 

- 

- 

He 

1.23 

1.00 

1.00 

0.417 

0.388 

0.397 

Li + 

0.763 

0.667 

0.667 

0.257 

0.248 

- 

Be ++ 

0. 569 

0.500 

0.500 

0.185 
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P (n £ I r) is taken to be the appropriate Hartree-Fock orbital. Dalgarno 
o i i 1 

and Parkinson (1959) have made the further modification of forming anti- 
symmetrical combinations of the product type solutions (81) and (90) 
and using them to calculate the polarizabilities. It is of interest to 
note that the application of this procedure to (114) with the variational 
solution (115) yields for the dipole polarizability 



( 120 ) 


a formula given by Heilman (1935), and application of it to the varia- 
tional solution (115) and the functional (114) yields 




( 121 ) 


a formula given by Buckingham (1937). 

We shall derive the correct version of (120) and (121) in Para. 
4.5. They differ because of the inconsistency involved in using a 
Hartree-Fock orbital in an equation derived from the Hartree approxima- 
tion scheme. This inconsistency is also reflected by the occurrence of 
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infinities in the potential derived from (113), wherever the Hartree- 
Fock orbital P (n.£. I r) has a node. In variational calculations, the 
difficulties presented by these unwanted singularities are suppressed 
by choosing trial wave functions of the form (119). 

In calculations using numerical integration, they are presumably 
treated by smoothing (r^) through the singularities (Dalgarno and 
Parkinson 1959). The presence of the singularities introduces some 
arbitrariness into the derived values of polarizabilities and shielding 
factors . 

4.4 THE STERNHEIMER PROCEDURE 

An instructive analysis of the Sternheimer procedure Is obtained 
by formally solving (94) as an expansion in terms of the eigenfunctions 
of (82), which we label more explicitly by the principal and azimuthal 
quantum numbers u^ = u (n£) (Sternheimer 1954). Then as in the deriva- 
tion of (52), we obtain 


a 


2L 


- Y. “21 < - nI) 


( 122 ) 


n Z 


where 


a 2L (ni) = -2 


3 1 (VlQ, v L u(nl)^ 


n ' 


e (o) (n'i') - e (o) (n Z) 


(123) 
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and 


2L 


' £ 7 2 L <"« 


(124) 


ni 


where 


7 2L ” ~ 2 




|u(n / i / ) , v T u(nJ)) (u(n'J'). v^ u(n £)j 
(n / £ / ) - (ni) 


(125) 


where the (ni) summation is over the occupied orbitals and the (n'i') 
summation is over all possible orbitals and includes an integration over 

the continuum. 

It is apparent from the form of (123) or (125) that in the evalua- 
tion of (106) - (109) the contribution to the polarizability or the 
shielding factor from a particular electron shell may be positive or 
negative but that in summing over all shells there will occur a nearly 
complete cancellation of contributions from transitions between occupied 
orbitals. The cancellation is not quite complete because (r ^ differs 
for different electron shells and u(ni) are not members of the same com- 
plete set. The Hartree approximation, with or without Sternheimer s 
modification, gives rise to spurious contributions from transitions 
between occupied orbitals. 
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Because of the energy denominators in (129) and (131) it is to 
be expected that upward transitions from the outermost shells will pro- 
vide the major contributions to the polarizabilities and shielding factors 
and this has been confirmed by explicit calculations (cf. Sternheimer 
1954, 1957 , Wikner and Das 1956). 

The summations over Z f in (123) or (125) can be written out expli- 
citly as in (103). Thus 


a d (ni) 




^ (n'i + 1), Vjl 
e (o) (n'i + 1) - 



(ni) 1 ) l 2 

(ni) 


= a. (ni - i + 1) + a (ni - i - 1), say (126) 

a q 


or 


(ni) 


-2 


u (0) (ni) v 2 u<°> (n'i + 2)J 


u (o) (n'i + 2), v 2 ' u (o) (ni)y 


e (o) (n'i + 2) - e (o) (ni) 


^ (ni), v 2 u (o) (n'i)^ (n'i), v 2 / u 

e ( °> (n'i) - e (o) (ni) 
n 

(127) 




=7 (ni - i + 2) + 7 (ni - i) + 7 ( ni - i - 2), say 

00 00 OO 
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As an illustration of the cancellation of contributions from 
transitions between occupied orbits, Sternheimer (1959) finds that for 
Na + , a d (2s - p) = 0.117A 3 and « d (2p - s) = -0.098A 3 , the largest contri- 
bution to a. (2s - P) and a (2p - s) arising from the 2s - 2p transition, 
da o3 

The most important contribution to a d is a d (2p - d) = 0.133A giving a 

o3 

net polarizability of 0.152A . 

For polarizabilities, transitions involving an increase in the 
azimuthal quantum number are usually the most important, especially for 
heavy ions with closed shell configurations. For Cs + (Sternheimer 1959), 
a (5p - f) = 5.21A 5 and a (5p - p) = 1.70A . For shielding factors, 

q q 

transitions involving no change in the azimuthal quantum number are 
usually much more important. For the quadrupole shielding factor of Na , 
Sternheimer (1959) obtains 0.30 for y m (2p - f) and -5.16 for 7<x> (2p - p) 
while for K + , he obtains -17.15 for y^Op ~ P)> -1-22 for ~ P) anci 

1.05 for the sum of the remaining terms all of which involve a change 

in i. 

For moderate and heavy ions then, the major contribution to y ^ 
comes from the radial modes of excitation y^ni - i)> the contribution 
from the angular modes yj ni -i+2) being comparatively small. The radial 
modes produce a negative shielding factor and so give rise to anti- 
shielding. A similar behavior occurs for hexadecapole moments (Stern- 
heimer 1961). 
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There is a real physical distinction between radial and angular 
modes of excitation which can be understood by considering a non-spherical 
nucleus placed at the center of a spherical shell of electric charge 
(cf. Kopfermann 1958). There will occur an angular redistribution which 
will concentrate the charge in the regions closest to the nuclear charge, 
screening the nuclear multipole moments and leading to a reduction in 
the total coupling. There will also occur a radial redistribution, the 
shell being distorted inwards by the electrostatic interaction between 
the positive nucleus and the negative electric charge. This radial 
redistribution gives an apparent increase in the nuclear coupling and 
is antishielding. 

4.5 THE UNCOUPLED HARTREE-FOCK APPROXIMATION 

The unphysical transitions between occupied orbitals can be excluded 
and the inconsistencies of the Sternheimer procedure can be removed by 
using total wave functions which satisfy the Pauli principle. To avoid 
mathematical complexity, we shall restrict further development to closed 
shell systems. 

The unperturbed wave function is written as 


x|r (r) 
T o ~ 



(128) 


46 



where 


4 


is the normalized anti symmetrizing operator and the unperturbed 

(o) 


orbitals u ^ ; satisfy the eigenvalue equations 


v t (0 V - h <0> “i (0) ^> 


(129) 


where 




= ( - 1/2 V 


(o). x I 2 

(. r k )l 


2 Z , r* r IV 

i - > u t fti> + L J TIT^ 

1 k=i 


dr. 


-f.J 


“ k (c> ^> “t gkj 4r u (o) 

, i K ^1 


k=l 


r . - r 


(130) 


the second summation being restricted to pairs of orbitals with parallel 
spins. The total unperturbed wave function ty Q (r) satisfies the equation 


(H' - E ') t (r) = 0 
N o o 


(131) 


where 


H / = 


z rv r | “, ( ° ) ^> > 2 k ( °V i 2 
L H i ' LL J I r , - r, I 


dr, dr 
~i ~k 


i<k 


zz J 

i<k 


- t (0> <£i> u k (0) ^> u i ( ° >( ^ ) u k <o)( ^> 


1 d ik 


r . - r. 

~i ~k 


(132) 
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and 



n „p 


(°) 

i v ~i 


(o) 


(o) 


(o) 


(?,) u^ vw/ (r 4 ) u/-(r k ) u k -^(r k ) 


dr, dr . 


i<k 


r . 


(133) 


Proceeding as for the uncoupled Hartree approximation, we solve 


(H 7 - E q O ^ (r) + v L f o (r) = 0 (134) 


by the substitution 


<£) ■ 



u. ^ (r. ) u. ^ (r . ) . 
i J ~J 


(135) 


Then (Dalgarno 1959) 
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~1 
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and 


« 2 L = -2 


Z f (u i <0)> v l u i (1>) ■ Z n <u 

i-l L k-l 11 


(o) - , . (o). 

k ■ V L ( Si ) “l > 


( Ul (1) , u k <»>, , 


(137) 


- 2 Z ( (u i <o) ' v £i> u i a)) - Z M <u k <o) - v4i> u i <o)) 

i=l V. k-l 11 


< u k 


( 1 ) 

> 



(138) 


If we now solve (136) formally as an expansion in the set of eigenfunc- 
tions of (131), we obtain expressions identical to (123) and (125) 
except that the ( n / £ / ) summation now excludes all the occupied orbitals. 

The solution can be achieved by methods similar to those 
employed in deriving the unperturbed orbitals, the substitution (98) 
again effecting a reduction to sets of radial equations. Alternatively 
a variational method can be employed based on the functionals 
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i v 1 ’ 1 1*1 


N 

V 

j=l 


(e. 
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The trial function 


(i) 


(r . ) = X v (r ) 


u.<°> (r ) 
1 ~1 


(140) 


yields the approximate formula for the dipole polarizability 
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where 


D. = 1 
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j=i 


(U J 


(o) 


(o). 




(°) , 


(°) 


(r, ) u ,, v ^ (?,• ) u t v^ k 
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first derived by Knipp (1939) from arguments based upon perturbation 
theory. The summation 


Sd 





(o). . (o). . (o) 





dr. 


d ^k 


(143) 


is the dipole oscillator strength sum rule appropriate to the eigenfunc- 
tions of the Hartree-Fock Hamiltonian (132). Knipp (1939) has also 
given a generalization of (141) which follows from (139) and the trial 
function (118), but it must considerably underestimate . 

The solution of (137) has been obtained only for beryllium 
(Dalgarno and McNamee 1961) and their results are shown in Table 3. 

The values obtained from the uncoupled approximation are much smaller 
than those obtained by the Sternheimer procedure and it seems that many 
of the shielding factors used in the interpretation of experimental data 
may be seriously in error. Beryllium is a special case, however, in that 
antishielding does not occur (within the Hartree-Fock approximation) 
whereas in cases of greatest experimental interest the antishielding 
is dominant. 
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TABLE 3 


2 2 

PROPERTIES OF B g (Is 2s ) S 



3 

V A) 

p oc (A) 

K oo q 

7 

' 00 

Sternheimer procedure 

9.5 

4.4 

15.1 

1.24 

Fully-coupled approximation 

4. 5 

1.8 

9.1 

0.77 


A modification of (136) has been used by Khubchandani , Sharma 
and Das (1962) to calculate the radial antishielding in C£ and K. In 
place of the exchange operator 




(144) 


k=l 


r. - r. 


they write, in effect, 


N r 

u 


u, (o) (r, ) u. (o) (r ) u (o) (r ) (l) 

JS ^ L dr — f ~ — u. (1) (r ). (145) 

-k (o) , „ r ~r 


k=l 


r . - r, 


u - v y (r.) 
1 ~1 
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The procedure can be regarded as an extension of the uncoupled 
Hartree approximation for determining u ^ ^ in which the inconsistencies 
arising from the use of unperturbed Hartree-Fock orbitals are eliminated. 
Khubchandani et al. find that compared to the predictions of the Stern- 
heimer procedure, the magnitude of 7 ^ for K is increased by 30% and for 
Cl' by 11% but there is an insignificant : -tcvatlon of the quadrupole 
polarizabilities . 

4.6 THE COUPLED HARTREE-FOCK APPROXIMATION 

The previous formulae (137) and (138) for polarizabilities and 
shielding factors remain valid in the fully-coupled Hartree-Fock approxi- 
mation, the only differences arising in the specification of the perturbed 
orbitals. The theoretical development has been given by Dalgarno (1959), 
Kaneko (1959) and Allen (1960). It follows from writing the wave func- 
tion t(r) of the perturbed system as 


N 

*<£> - n 


and minimizing the total energy with respect to arbitrary variations of 

u ^ and u Then expanding in powers of the perturbation, the first 

i i 

order perturbed orbitals satisfy the equation 


u. (o) (r.) + u. ( 1 ) (r.) + 
i i ~i 


(146) 
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( H i - £i (0) ) u. (1) (r.) + v L (r.) - 6^“' ju. 


U. (0) (r.) 


' - Vws'-V ■•>„., 

I ~k i ~1 

r, - r 

1 


N r 

k=l 


N 

z 

k=l 



Js 2* i 2*. dr, U, ( 1 ) (r.) + 

, ~k k ~r 

r. - r, 

~i ~k 1 
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dr, a, (o) (r. 
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(147) 


and 


N (!),_ x .. (o),_ \ | .. (o), 

^ (1> - <» 1 <0> . V£ t > + 2 I / 
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In (147), the perturbed orbitals are coupled together both by 
direct and exchange interactions. The substitution (98) reduces (147) 
to coupled sets of radial equations but the solution involves laborious 
computation and it has been carried out only for the helium sequence. 

The fully-coupled approximation is, however, much more accurate than 
the uncoupled approximation. Thus, if the difference between the Hartree 
Fock Hamiltonian (132) and the actual Hamiltonian of the system is 
regarded as a small perturbation, then xn the absence of degeneracies 
the polarizabilities derived from the fully-coupled approximation are 
correct to first order whereas for the polarizabilities derived from 
the uncoupled approximation there is a non-vanishing first order correc- 
tion (Dalgarno 1959). By an obvious extension of the arguments given 
by Cohen and Dalgarno (1961), it may be shown that the theorem applies 
also to the calculation of shielding factors. In addition, the fully 
coupled approximation will usually yield lower bounds for the polarizabil 
ities (Dalgarno and McNamee 1961a). 

For the helium sequence, (147) simplifies to 


( - 1/2 V. 2 + V.(r.) - e. (0) ) u (1) (r.) + v L (r.) 


/ 


•‘"fa* 

r, - r 

1 


(149) 
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Solutions of (149) for dipole and quadrupole perturbations have 
been obtained by Dalgarno (1959), Kaneko (1959) and Dalgarno and McNamee 
(1961) and the resulting values of the polarizabilities and shielding 
factors are given in Table 2. The results are encouraging (again with 
the exception of H~), especially for the quadrupole shielding factor. 

The dipole shielding factors of the helium sequence are actually given 
exactly by the fully-coupled approximation and serve merely as a comput- 
ational check. 

Since the effect of correlation is usually to reduce the polariza- 
bilities, the uncoupled Hartree approximation probably overstimates 

as it does a,. Accordingly, for He 
d 

o5 o5 

0.0965A < a < 0.0979A (150) 

- q - 


and for Li + 


o5 o5 

0.00464A < a < 0.0047 2A . (151) 

- q - 

The lower bound for helium is larger than the value given by the more 
refined variational calculation of Dalgarno, Davison and Stewart (I960), 
suggesting that their adopted six-parameter representation of the unper- 
turbed wave function does not provide an accurate description of the 
electron density at large distances from the nucleus. 
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Watson and Freeman (1962) have pointed out that calculations of 
shielding factors can be carried out within the framework of the Hartree- 
Fock method, provided two constraints are relaxed. Relaxing the con- 
straint that the orbitals be separable into products of radial and angular 
factors gives the angular excitation contribution and relaxing the con- 
straint that the radial wave functions be the same for electrons with 
different magnetic quantum numbers in a particular electron shell gives 
the radial excitation contributions. They have applied their method to 
the calculation of the radial excitation contributions to the quadrupole 
shielding factors of Cl", Cu + , Br", Rb + , Cs + , l“ and Ce 3+ . For Cu + , 
they obtain - 17.6 in good agreement with the value of - 17.1 computed 
by Sternheimer and Foley (1956) but for Cl", they obtain - 88.9 which is 
much larger in magnitude than the value of - 56.5 computed by Sternheimer 
and Foley.* It would be interesting to compare the results of Watson 
and Freeman with those given by the fully-coupled approximation. 


X am indebted to Dr. A. Freeman for providing me with this description 
of his work in advance of publication. 
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SECTION 5 


EXPANSION METHODS 

An assessment of the accuracy of the various methods of calculating 
polarizabilities and shielding factors can be obtained by expanding the 
wave functions as series in inverse powers of the nuclear charge, a pro- 
cedure which uncouples the perturbed orbitals and leads to equations 
which can be solved analytically. Its usefulness in predicting polariz- 
abilities and shielding factors depends upon the rapidity of convergence 
of the derived series expansions. In general, the convergence will be 
more rapid for shielding factors, which depend upon the distortion near 
the nucleus, than for polarizabilities, which depend upon the distortion 
far from the nucleus. The method has the advantage that results for all 
members of an isoelectronic sequence are obtained in a simple calculation. 

5.1 THE HELIUM SEQUENCE 

Expansion methods have been used by Dalgarno and Stewart (1958) and 
Dalgarno , Davison and Stewart (1960) to solve the exact Schrodinger 
equations appropriate to the dipole and quadrupole polarizations of the 
helium sequence. They expand 

* - E ♦ (S) , E = E E (S) (152) 

T n s n n s n 
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where 


4 > 
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<£i 


~ 2 } 


(Z /it) exp 




( 153 ) 
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o 



( 154 ) 


the unperturbed Hamiltonian being 


H" = - 


1 

2 



Iv 2 
2 2 




( 155 ) 


To first order in Z , 
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2L 


= 2 




) 



( 156 ) 


(o) 

with a similar expression for 7 2 L * The first order functions <t>^ , 

d> ^ and 4> ° are solutions of 
o o 


(H" 




= 0 


(157) 
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(H // 



\ 



( 158 ) 







* ° = 0 . ( 159 ) 

o 

The solution of ( 158 ) can be obtained immediately but, because 

of the electron interaction term, the solutions of ( 157 ) and ( 159 ) pre- 
sent great difficulty. However, using the methods of Dalgarno and Stewart 
( 1956 ), ( 156 ) can be written in the alternative form 
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where 


(H 


/ / 
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(°) 


+ < V L 



(o) 


- (Cl v l -Ei (o) |C o) ) C = ° (161) 

and the solution of (161) can be achieved by elementary means. The 
analogous expression for 72 ^ is (Dalgarno, Davison and Stewart 1960) 


7 2L 
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(162) 


where 


CH" - E o ) 


'11 + <V, L - E 1 <0),) *1 <0) + (U 1 - E 1 (0) 


) * 


(o) / 




(163) 
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and the solution of ( 163 ) can also be achieved by elementary means. 

The details can be found in the papers cited. The results are that 
to first order in Z 


a 


d 



2L, + o(z- 6 ) 

16Z 


(164) 


2 

Z 


(165) 


a 

q 



+ 0(z - 8) 

128z' 


(166) 


7 = + — - — ~ (480 in 2 - 329) + 0(z" 3 ) (167) 

3Z 18Z 2 

Applying the screening approximation which takes some account of higher 
order terms (Dalgarno and Stewart 1960), we obtain 


a. = 9/(Z-0.3594 ) 4 (168) 

d 

B = 2/Z (169) 

'CO 

a = 30/(Z-0.4506) 6 (170) 

q 

7 = 2/3 (Z-0. 3092) . (171) 

* CO 
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The accuracy of the simple formulae is superior to that of the 
uncoupled Hartree approximation and is comparable to that of the coupled 
Hartree approximation. Thus (171) gives a value of 0.394 for the quad- 
ruple shielding factor of helium compared to the (probably) accurate 
value of 0.396. The accuracy of the formulae increases rapidly with 
increasing nuclear charge. 

A similar analysis can be applied to the Hartree-Fock approximation 
(Cohen and Dalgarno 1961). For the helium sequence, it is found that 
the uncoupled approximation is equivalent to the formulae 


a 


d 



1923 

128Z 5 


o(z" 6 ) 


P 


00 



0(z" 3 ) 


(Cohen and Dalgarno 1961, unpublished) 


a 

q 


30 10485 

Z 6 128Z 7 


0(z" 8 ) 


3Z 


13 

48Z' 


+ 0(z“ J ) 


(172) 

(173) 

(174) 

(175) 


(Cohen, Dalgarno and McNamee 1962), whereas the coupled approximation is 
equivalent to formulae identical to (168)", (169), (170) and (171). The 
coupled approximation therefore leads to results which are correct to 
first order in harmony with the assertions in Para. 4.6. For (3^ the 
coupled approximation actually leads to results which are correct to all 

orders of Z 


63 


5.2 THE BERYLLIUM SEQUENCE 


Cohen, Dalgarno and McNamee (1962) have also conducted an analysis 
of the calculation of the quadrupole polariEabilities and nuclear shield- 
ing factors of the beryllium sequence using the uncoupled and coupled 
Hartree-Fock approximations. According to the uncoupled approximation, 
they find that 


a 

q 


7 

' 00 


32670 


(1 


8 . 002 . 

7 ' 


2 

Z 


(1 




(176) 

(177) 


and according to the coupled approximation they find that 


a 

q 


32670 


9 199 

(i + ^|^) 


1 (i + i#2) 


(178) 

(179) 


It appears from these comparisons that only a small improvement in 
accuracy will result from using the coupled approximation in place of 
the uncoupled approximation. It will be of interest to extend the com- 
parison to cases where antishielding can occur. 

The series expansions for ct converge too slowly to be quantitatively 
useful but we may obtaip values of y ^ of acceptable accuracy from (177) 
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and (179). Applying the screening approximation (Dalgarno and Stewart 
I960), we obtain for the uncoupled approximation 

7 = 2/ (Z- 1 • 263) (180) 

which yields a value of 0.73 for beryllium in good agreement with the 
value of 0.77 obtained by numerical integration by Dalgarno and McNamee 
(1961b). The coupled approximation yields the more accurate formulae 

7 = 2/ (Z-l . 220) (181) 

' CO 

according to which y^ is 0,72 for beryllium. 

Contrary to the case of the helium sequence, (178) and (179) are 
not identical to the exact expansions. This occurs because the Hartree- 
Fock approximation does not take account of the degeneracy in the limit 
of infinite Z of the substates common to a particular quantum number 
(Layzer 1959, Linderberg and Shull 1960). Thus the correct zero order 
wave function for the ground states of the beryllium sequence is not 
\|/ (Is 2 2s 2 ) 1 S as the Hartree-Fock approximation asserts but Jne inear 

combination 


(^S) = a ^(ls 2 2s 2 ) 3 + b ^ Q (ls 2p ) S (182) 
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where a and b are constants so chosen that the matrix of the Hamiltonian 
is diagonal. The modification is important for the second component of 
the wave function introduces the possibility of antishielding. Cohen 
et al. (1962) have calculated the zero order changes in the quadrupole 
polarizability and shielding factor and they find that the leading term 
of (178) should be 31541/Z b in place of 32670/Z and that the leading 
term of (179) should be 1.784/Z in place of 2/Z. It appears that the 
coupled approximation may overestimate o:^ by about 4 L and by about 
11 %. 

In the case of the dipole polarizability of beryllium the calcula- 
tion would be further complicated because the dipole perturbation itself 

2 1 

couples the two components of (182) to ^(ls 2s2p) P, with which they 
are degenerate. To obtain reliable values of a d> it seems necessary to 
adopt a multi-configuration approach of the kind explored by Donath (1961). 
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SECTION 6 


STATISTICAL METHODS 

Calculations of polarizabilities (Gombas 1944, 1956) and of shield- 
ing factors (Sternheiraer 1950, 1951, 1954, 1957, 1961; Sternheimer and 
Foley 1953) have been carried out based upon the Thomas -Fermi model of 
the atom. A comparison of the dipole polarizabilities computed by Gombas 
(1944) with the most accurate values available (cf. Table 5) is given 
in Table 4. The method overestimates a d considerably for neutral atoms 
but its accuracy improves with increasing mass and with increasing excess 
charge . 

A similar behaviour may be anticipated for shielding factors but as 
currently employed the statistical description does not yield any informa- 
tion on the radial antishielding which is dominant for heavy systems. 
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TABLE 4 


-24 3 

DIPOLE POLARIZABILITIES IN UNITS OF 10 cm 


System Statistical Theory Accurate Value 


f" 

6.20 

-1.2 

Cl" 

7.10 

-3 

Br~ 

8.41 

-4.5 

i’ 

9.21 

-7 

Ne 

2.01 

0.40 

A 

2.88 

1.64 

Kr 

4.00 

2.48 

Xe 

4.61 

4.04 

Ne + 

0.850 

0.15 

K + 

1.36 

-0.9 

Rb + 

2.14 

-1.7 

Cs + 

2.66 

-2.5 

Mg ++ 

0.400 

0.07 

Ca^ 

0.721 

-0.5 

Sr^ 

1.30 

-1.1 

Ba" 1- *" 

1.70 

-1.7 
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TABLE 5 W 


DIPOLE POLARIZABILITIES IN UNITS OF 10 


System 

(a) 

(b) 

System 

(a) 

(b) 

H~ 


31.4< 3 > 

C 

2.1< 21 > 



13.4< 2 > 

30.2 (4) 

o 2+ 

0.3< 16) 


He 

0.218 (1) 

0. 196 ( ' 5,6 ^ 

Ne 3+ 

0.13< 16) 



0. 224^ 2 ^ 

0. 2050 ( - 3 ' ) 

N 

1.3< 21) 

i.i 3 < 22) 


0.220^ 

0.2051 (?) 

0 + 

0.49< 16) 


Li + 

0.0305 (1) 

o.0286^ 8) 

0 

0.89^ 21 ^ 

0.77< 22 > 


0.0307 (2) 

0.0281 (5) 

Ne 2+ 

o.i 5 < 16) 



0.0304 (5) 

0.08^ 9,10 ^ 

0“ 

3.2< 16 > 




0.024< U ’ 12) 

F 

o.6< 21 > 


Be 2+ 

0.00813 ( ' 1) 

0.00759 (14) 

Ne + 

0.21^ 



0.00825 (2) 

0.04< 9) 

f" 

1.9< 23) 

1.0< 9 ' 10) 



0.007( U > 12 ) 


1.8< 21 > 

0.76< 13 > 

b 3+ 

0.00303 (2) 

0.00288 (14) 



1.2< 26 > 



0.02< 9 > 

Ne 

0.40 (21) 

0.395^ 7 ' 



0.003^ U > 12) 


0. 62^ 6 ^ 

0.367 <26 ^ 

c 4+ 

0.00139 (2) 

0.00132 (14:) 

Na + 

0.152< 23 > 

0.2< 9 - 10) 



0.0013 (8) 


0.145 (16 ' ) 

0.26< 13 > 



0.012 (9) 



0. 154* 26) 
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(9,U) 


Li 

2C< 15) 

24.4< 1S > 

Mg 2+ 

0.082 (L6) 

0.10< 9 ’ U > 


25 (16) 

20 ± 3 (19) 



0.072 ( ' 24) 


24.9*- 17 -* 

16 + 3 ( ^ 25 ' ) 



0.012< 10 > 

Be 

9.6< 16 > 

4. 5 ( 20) 

ai 3+ 

0. 050^ 12 ^ 

0.053 (11) 


9.5< 20) 



0.045< 16 > 

0 . 07 ^ 9 5 10 ^ 

B 

5.1 (21) 


Si 4+ 

0.027 (16 ' ) 

0.04< 9 ’ 10 > 

o 3+ 

0.3< 16) 


Na 

22.9 (1?) 

24. 6 






20 + 2.5 (19) 

ci" 

7.i9 (23) . 

3.5< 9) 

2+ 

Sr 


1.4< 10 > 



3.0^ 10 ’ 13) 



1.6< L3 > 

A 

2.40< 6 > 

1.64< 7) 

Rb 

49.1< 17) 

43.9< 18 > 

K + 

1.24< 23 > 

0 <9 (9,10) 



40 ± 5< 19 > 



o.8< u > 

^ 2+ 
Ca 


1.8< 13) 



1 . 2 ( 13 ^ 

_ 4+ 
Sn 


3.4< 13 > 

Ca 2+ 

0.73 (12 ' 15) 

0.54< U > 9 > 

i" 


7 . 6 ^ 9 ^ 

„ 3+ 
Sc 


0.35^ 9) 



6.2< 10 ’ 13) 

^.4+ 

Ti 


0. 24^ 

Xe 


4.04 (7) 

K 

44.4< 1? > 

41.6< 18 > 

Cs + 

5.60< 23 > + 

2.7 (9 ’ 10) 



40 + 5< 19) 



2.4( U ) 

Cs 

57 ( 15 ) 




3.1 (13) 

Sc 

154< 15) 


Ba 2+ 


1.7< 9) 

Fe" H ' 


_ 1 (13) 



2.5< 13) 

Cu + 

0.982 (23 ^ 

i.6< 13 > 

T 3 + 

La 


1.3< 9) 

Zn 2+ 


0.8< 13 > 

Cs 

67.7< 17) 

53.8^ 
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„ 4+ 
Ge 


i.o (13) 



52.5 +6.5 

2+ 

Cu 


o.2< 13 > 

Ba 


62< 25 > 

Br' 


4.97 (9) 

pb 2+ 


4.9< 13 > 



4.15 (10)13) 

Hg 2+ 

2.78 (23)t 


Kr 


2.48 (7) 

Te + 


5.2< 13 > 

Rb + 

2.92< 23)+ 

1 . 7 (9 ’ 10) 

u 6+ 

1.34< 23 * 




1.8< 13) 




Sr 2+ 


0.86 (9) 






1.0< U > 





Footnotes for Table 5 


^Values for a few systems not included in Table 4 have been estimated by 
Pauling (1927) using a screened hydrogenic formula. 


+ These values are computed using Hartree wave functions, which usually 

overestimate CL , considerably, 
d 

(a) Values computed using the Sternheimer procedure. In References (16) 
and (21), the antisymmetrical version has been used. 

(b) Experimental values and more accurate theoretical values. 
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SECTION 7 


NUMERICAL VALUES 


7.1 DIPOLE POLARIZABILITIES 

A compilation of dipole polarizabilities computed using the Stern- 
heimer procedure or the antisymmetrical version of it is given in Table 
5. When the result of a more refined theoretical calculation exists it 
also is included. 

Experimental values of dipole polarizabilities can be obtained from 
deflection measurements such as those of Salop, Pollack and Bederson 
(1961), from extrapolation of measurements of refractive indices (cf. 
Dalgarno and Kingston 1960), from analysis of the Rydberg-Ritz correc- 
tions for spectral series (Born and Heisenberg 1924; Mayer and Mayer 
1933), from the indices of refraction of salts in aqueous solutions 
(Fajans and Joos 1924) and from refraction data on crystals (Tessman, 
Kahn and Shockley 1953). Table 4 includes the experimental values. 

An assessment of the accuracy of the theoretical values is confused 
by defects in the theories used to interpret the experimental data and 
by environmental effects which cause the polarizability to differ from 
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that of the free atom or ion (Burns and Wikner 1961). The Sternheimer 
procedure, being essentially a single-electron approximation, should 
yield accurate dipole polarizabilities for the alkali metals but its 
accuracy is unlikely to be sustained for many-electron systems, especially 
those with an outermost shell of electrons described by orbitals with 
nodes other than at the nucleus. Table 5 is, at least, not in disagree- 
ment with this view. It appear that with the exceptions of beryllium 
and boron the polarizabilities predicted for neutral atoms lighter than 
sodium are unlikely to be in error by more than 15%, a figure also 
suggested by some calculations of the induced electric field at the 
nucleus by Sternheimer (1959). For beryllium, boron and atoms heavier 
than sodium, the possible error may approach a factor of two. For 
positive ions the errors are smaller and for negative ions they are 
larger. It is usually argued that correlation effects act to reduce 
polarizabilities, suggesting that the predicted polarizabilities are 
too large. This is probably true in general but not always, as the 
example of H demonstrates. 

The table provides clear evidence that environmental effects signi- 
ficantly modify polarizabilities. Thus the dipole polarizability of 

, -24 3 

free Na^ cannot be very different from 0.15 x 10 cm whereas the 

-24 

value derived from the refractivity of ionic crystals is 0.26 x 10 
3 

cm . 
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7.2 QUADRUPOLE POLARIZABILITIES 


A compilation of quadrupole polarizabilities computed using the 
Sternheimer procedure is given in Table 6. The accuracy of the predicted 
quadrupole polarizabilities is comparable to but somewhat higher than 
that of the predicted dipole polarizabilities. 

7.3 QUADRUPOLE SHIELDING FACTORS 

Quadrupole shielding factors, computed according to the Sternheimer 
procedure, are collected together in Table 7. The results of more 
refined theoretical calculations are also included and are given in 
brackets. The accuracy of the predicted values is comparable to but 
somewhat lower than that of the dipole polarizabilities. 

The significant feature of the table is the large antishielding 
predicted for heavy systems, to which Foley, Sternheimer and Tycko (1954) 
first drew attention. Qualitative confirmation of the existence of 
large antishielding effects has been obtained from nuclear quadrupole 
coupling data in alkali-halide gases (cf. Townes 1958; Wikner and Das 
1958; Burns 1959a) and in ionic solids (cf. Bersohn 1958; Bernheim and 
Gutowsky 1960; Simmons and Slichter 1961; Hewitt and Taylor 1962), from 
analysis of nuclear spin transitions in crystals induced by ultrasonic 
waves (Proctor and Robinson 1956; Jennings, Tanttila and Kraus 1958; 

Taylor and Bloemberger 1959; Bolef and Mernes 1959), from nuclear magnetic 
relaxation times in ionic crystals (Van Kranendonk 1954; Das, Roy and 
Ghosh Roy 1956; Wikner, Blumberg and Hahn I960), from dislocation studies 
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TABLE 6* 

* 

-40 5 

QUADRUPOLE POLARIZABILITIES IN UNITS OF 10 cm 


System 


System 


h" 

66.5< 27 > 

Mg 2+ 

0.0223 (30) 


67.0< 28) 

ai 3+ 

o.oioi (28) 

He 

0.0993 (27) 


0.00915 (30) 


0.0949 (28) 

Si 4+ 

0.00438 (30) 


0.0979 (5) 

Cl" 

13.8 (27:> 

Li + 

0.00473 (27) 


13.1 (29)+ 


0.00471 (28) 

A 

2.19< 30) 


0.00472 (5) 

K + 

0.733 (27) 

Be 2+ 

0.000637 (27) 

Ca 2+ 

0.721< 31 >* 


0.000633 (28) 

Cu + 

0.0309 (30) 

b 3+ 

0. 000142 ( ' 28 ' ) 

Rb + 

1.28 (27) 

Be 

15.1< 20) 


2.99 (27)+ 

f" 

2.38< 29 > 

Cs + 

3.03< 31)+ 

Ne 

0.370< 3 °> 


7.80< 27 > + 

Na + 

0.0670 (28 ' ) 


7.86< 31 > t 


0.0634 (23) 




0.0649 (31) 



* i 

A partial 

evaluation of CL has been 

carried out 

by Wikner and Das (1958) 

for Cl , 

K + , Cu + , Ga 3+ , Rb + , Cs + , Br" and i". 


^ Computed 

using unperturbed Hartree 

orbitals 



* These values are in harmony with the results of a more refined calcula- 
tion by Khubchandani e t al . (1962). 
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TABLE 7 


QUADRUPOLE SHIELDING FACTORS 1 " 


h" 

i.i«< 28) 

f" 

-25.7< 29 >* 

Ca 2+ 

- 8.4 (30)t 


i.m< 23 > 

Ne 

- 8.99 (30) * 

Mn 2+ 

- 6 . 81 ^ 29 ^ 

He 

0. 416 ( ' 28 ^ 

Na + 

- 4.1 (34 > 

Fe 3+ 

- 6.17< 29 > 


0.424 (23;> 


- 4.5< 28) 

Cu + 

-15.0< 34 > 


(0.397) (32) 


- 4.6< 23 > 


-13.8 (35) 

Li + 

0.056 (28) 

Mg 2+ 

- 2.7< 30 >* 


(-15.1) (36) * 


0. 263^ 23 ' ) 

ai 3+ 

- 2.6< 28 > 

Ga 3+ 

- 6.94< 29 > 


(0.249) (32) 

si 4+ 

- 1.9< 30) * 

Br" 

-97.0 (35 >** 

Be ++ 

0.185< 28) 

Cl" 

-56.6^ 

Kr 

-41.0< 29) 


0.189< 23) 


-50.1< 29) 

Rb + 

-49 . 3 O 5 )** 


(0. 181) (32) 


-49.3< 35) 


(34)** 
-70. 7 V ' 

b 3+ 

0.145< 28) 


(-55.5) (3?) 

cd 2+ 

_ 15 (29)tt 


(0.142) (32) 


(-87.5)< 36) * 

Ag + 

-22.1< 29) 

Be 

1 . 24^ 20 ^ 

A 

-21.2< 30) * 

T 3 + 

In 

-15.3 (29) ** # 


(0.72) (33) 

K + 

-17.3< 23) 

i” 

_179 ( 3 5) ** 

B + 

0.768< 23) 


-12.8< 35) 

Cs + 

( 35 )** 
-11CT ' 


(0.53) (33 ' > 


(-16.3) (37) 


-144^ ; 

t A 

partial evaluation 

of for 

Cu + , V++, Ag + 

and Hg"*""^ has 

been reported 

rece 

ntly by Sternheimer 

(1962b). 




* The outer shell contribution on 

ly has been included. 



* In (36), the radial contribution only has been calculated. The value 
given includes the angular contribution obtained in reference (34). 

** These values are computed using Hartree wave functions. 

tt This is an interpolated value. 

# In reference (29), it is suggested that the use of Hartree-Fock wave 
functions for In 3 + would give a value of -11. 
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References to Tables 5, 6, and 7. 


(1) Wikner and Das (1957) 

(2) Sternheimer (1957) 

(3) Schwartz (1961) 

(4) Geltman (1962) 

(5) Dalgarno and McNamee (1961a) 

(6) Kaneko (1959) 

(7) Dalgarno and Kingston (1960) 

(8) Baber and Hasse 7 (1937) 

(9) Fajan and Joos (1924) 

(..-?) Born and Heisenberg (1924) 

(11) Mayer and Mayer (1933) 

(12) Sternheimer (1954) 

(13) Tessman, Kahn and Shockley (1953) 

(14) Dalgarno and Stewart (1958) 

(15) Sundbom (1958) 

(16) Parkinson (1960) 

(17) Sternheimer (1962) 

(18) Dalgarno and Kingston (1959) 

(19) Salop, Pollack and Bederson (1961) 

(20) Dalgarno and McNamee (1961b) 

(21) Dalgarno and Parkinson (1959) 

(22) Alpher and White (1959) 
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References to Tables 5, 6, and 7. (Continued) 


(23) 

S ternhe imer (1959) 


(24) 

Bockasten (1956) 


(25) 

Drechsler and Muller (1952) 


(26) 

Donath (1961) 


(27) 

S ternhe imer (1957 ) 


(28) 

Das and Bersohn (1956) 


(29) 

Burns and Wikner (1961) 


(30) 

Burns (1959b) 


(31) 

Burns (1959a) 


(32) 

Dalgarno, Davison and Stewart (1960) 

(33) 

Cohen, Dalgarno and McNamee 

(1962) 

(34) 

Sternheimer and Foley (1956) 


(35) 

Wikner and Das (1958) 


(36) 

Watson and Freeman (1962) 


(37) 

Khubchandani, Sharma and Das 

(1962). 
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(Otsuka 1958), and from mixed crystal experiments (Kawamura, Otsuka and 
Ishiwatari 1956; Otsuka and Kawamura 1957). There are serious deficien- 
cies in the theories used to interpret the experiments. Some measure of 
quantitative agreement is found between the predicted and derived values 
of y for positive ions but the predicted values of 7 for negative ions 
are apparently too large by a factor of about four (Burns and Wikner 1961). 
Although errors of this order cannot be excluded, the work of Watson and 
Freeman (1962) suggests that the Sternheimer procedure underestimates 
I y I and modification of the negative ion structure by its environment 
must be a major source of the discrepancies (Burns and Wikner 1961). 

7.4 HIGHER ORDER POLARIZABILITIES AND SHIELDING FACTORS 

Some calculations of higher order polarizabilities have been 
reported by Burns (1959a) and by Stewart (1961) and of higher order 
shielding factors by Sternheimer (1961). The hexadecapole antishielding 
factors are very large and higher order effects are important (Stern- 
heimer 1962b). 
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